Superconducting Gap Structure of LaFePO Studied by Thermal Conductivity 
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The superconducting gap structure of LaFePO (T c = 7.4 K) is studied by thermal conductivity 
(/t) at low temperatures in fields H parallel and perpendicular to the c axis. A clear two-step field 
dependence of k{H) with a characteristic field H s (~ 350 Oe) much lower than the upper critical 
field H C 2 is observed. In spite of large anisotropy of H C 2, m(H) in both if-directions is nearly 
identical below H 3 . Above H 3 , k{H) grows gradually with H with a convex curvature, followed 
by a steep increase with strong upward curvature near H C 2- These results indicate the multigap 
superconductivity with active two-dimensional (2D) and passive 3D bands having contrasting gap 
values. Together with the recent penetration depth results, we suggest that the 2D bands consist of 
nodal and nodeless ones, consistent with the extended s-wave symmetry. 

PACS numbers: 74.25.Fy, 74.20.Rp, 74.25. Op, 74.70.-b 



Recent discovery of a new class of Fe-based supercon- 
ductors pj has attracted much attention. Among them, 
FeAs-based compounds have aroused great interest be- 
cause of the high transition temperature T c . Undoped 
arsenide LaFeAsO(La-llll) is non-superconducting and 
has a spin-density-wave (SDW) ground state, but be- 
comes superconducting (T c — 25 K) when electrons are 
doped Q. By changing the rare-earth ion, T c reaches 
as high as 55 K in SmFeAs(0,F) Q. A key question is 
the origin of the pairing interaction. Since the symmetry 
of the superconducting order parameter is intimately re- 
lated to the pairing interaction at the microscopic level, 
its identification is of primary importance. 

Fully gapped superconducting states in FeAs-based su- 
perconductors have been reported by the penetration 
depth measurements of Pr-1111 Q, Sm-1111 [5], and 
Bai_ K K x Fe2As2 (Ba-122) [6|, angle-resolved photoemis- 
sion |3, thermal conductivity 0] and NMR relaxation 
rate [9|] measurements of Ba-122. Some of them give ev- 
idence of multiband superconductivity with two distinct 
gaps. On the other hand, the NMR of La-1111 [nj and 
Pr-1111 [ll[ and the penetration depth measurements of 
Ba(Fei_ 2: Co 2: )2As2 [12| suggest the presence of low-lying 
excitations, which could be indicative of nodes. Theo- 
retically, it is proposed that a good nesting between hole 
and electron pockets prefers the "s±" symmetry where 
the gap is finite at all Fermi surfaces but changes its sign 
on different bands [13j |. Recent neutron resonant scat- 
tering [14| and the impurity effects on the penetration 
depth [y, are consistent with this symmetry. 

The phosphide LaFePO isomorphic to LaFeAsO ex- 
hibits superconductivity at T c ~ 7 K without doping [lj| . 
In spite of the similarity in the overall electronic struc- 
ture [17| , the magnetic and superconducting properties of 
LaFePO and LaFeAsO are very different. In fact, the un- 
doped stoichiometric LaFePO is non-magnetic. Recently, 
a superconducting state with well-developed line nodes in 
the gap function has been suggested in LaFePO by a lin- 
ear temperature dependence of the penetration depth at 



low temperatures [LSI ]. However, on which Fermi surfaces 
the nodes locate in the multiband electronic structure is 
not yet clarified, and while several candidates have been 
theoretically proposed [H, H3| , the superconducting sym- 
metry in LaFePO remains elusive. Thus the clarification 
of the detailed gap structure of LaFePO is expected to 
provide important clues to the origins of magnetism and 
superconductivity of Fe-based compounds. 

Here, to shed further light on the gap symmetry of 
LaFePO, we present the thermal conductivity measure- 
ments at low temperatures. The thermal conductivity 
probes delocalized low-energy quasiparticle excitations 
and is an extremely sensitive probe of the anisotropy of 
the gap amplitude. We provide strong evidence of the 
multigap superconductivity, in a more dramatic fashion 
than FeAs-based superconductors, with two very differ- 
ent gap values. We show that there are passive 3D bands 
and two kinds of active 2D bands; one is fully gapped 
and the nodes inferred from the penetration depth mea- 
surements [l8| are most likely on the other 2D bands. 
This is compatible with the extended s-wave (nodal s±) 
symmetry for the gap structure of LaFePO. 

Single crystals with typical dimensions of 0.8 x 0.4 x 
0.05 mm 3 were grown by a Sn-flux method [21| . We care- 
fully removed the Sn-flux at the surface of the crystals 
by the diluted hydrochloric acid. The resistivity and sus- 
ceptibility measurements show the sharp superconduct- 
ing transition at T c = 7.4 K (determined by the midpoint 
of the resistive transition) , which is slightly higher than 
the values reported by other groups [lg, [22| • The ther- 
mal conductivity n was measured by a standard four- wire 
steady method for a heat current q within the ab plane. 

The temperature dependence of the in-plane resistiv- 
ity p in the zero field (inset of Fig. 1) depends on T 
as p = po + AT 2 with po — 4.9 /if2 cm and A — 
3.3 x 10~ 3 /jficm/K 2 below 50 K down to T c . The 
residual resistivity ratio (RRR) is 28. Clear de Haas- 
van Alphen (dHvA) oscillations were observed in sam- 
ples from the same batch with nearly the 
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FIG. 1: (color online) Temperature dependence of k/T in 
zero field and at poH = 9 T in the normal state above H^2 
(if || ab, HA. q). Inset: p plotted as a function of T 2 . 



value [2l| . The upper critical fields at T — > K estimated 
by the resistivity measurements are p,oH^ 2 = 1.0 T for 
H\\ c axis and n H^ = 8.6 T for ff|| ab plane. This 
large anisotropy of the upper critical field H C 2 indicates 
that the bands active for the superconductivity have a 
very anisotropic 2D electronic structure. 

Figure Q] depicts the temperature dependence of k/T 
in zero field and in the normal state above Hjg. As 
the temperature is lowered, k/T decreases below T c . 
The electronic contribution in k at 1 K estimated by 
the Wiedemann-Franz law, k = LqT / p (Lq — 2.44 x 
1(T 8 flW/K is the Sommerfeld value), is 0.5 W/K 2 m, 
which is close to the observed value 0.65 W/K 2 m. Thus 
the electron contribution dominates well the phonon heat 
contribution at least below 1 K. 

First we discuss the thermal conductivity in zero field. 
A residual term kqq/T at T — > K in k/T is clearly 
resolved. In the nodal superconductors, such a resid- 
ual term appears as a result of the impurity scattering 
which induces quasiparticles even at T — K. In the 
presence of line nodes in a single band superconductor, 
kqo/T is roughly estimated as ~ 2(^ a i,/£) ■ (k„/T), where 
£ a b is the in-plane coherence length, I is the mean free 
path and K n is the thermal conductivity in the normal 
state. Using I = 94 nm from the dHvA measurements 
[HI and i ab = ^$0/(2^3) = 18 nm, K m /T is esti- 
mated to be ~0.19 W/K 2 m. This value is comparable 
to the observed kqq/T ~ 0.30 W/K 2 m, but we note that 
this estimate includes large ambiguities due to the multi- 
band effect which could alter effective £, a b and K n . So 
this comparison alone cannot be taken as the evidence for 
line nodes in the superconducting gap. It should be also 
noted that the residual term may arise from an extrinsic 
origin, such as non-superconducting metallic region [22| 
with high thermal conductivity, though the sharp super- 
conducting transition and the observation of the dHvA 
oscillation indicate a good quality of the crystal. 



More vital information on the gap structure can be 
provided by field dependence of thermal conductivity at 
low temperatures. Field dependent part of k{H) in a 
mixed state mainly stems from the superconducting part 
of the crystals, even if a non-superconducting region was 
present in the crystal. Moreover, the phonon scattering 
at the low temperatures is governed by static defects and 
is therefore field independent. 

Figures 2(a) and (b) depict the field dependence of 
k(H) — k(0) normalized by n(H C 2) — k(0) for H\\ c and 
H\\ ab, respectively, measured at T = 0.46 K (0.062T C ) 
by sweeping H after zero-field cooling. We note that lit- 
tle difference was observed between the data measured in 
zero-field and field cooling conditions, indicating that the 
field trapping effect is very small. For both field direc- 
tions, the overall ^/-dependence of thermal conductivity 
is quite similar. At very low fields, the thermal conduc- 
tivity exhibits a pronounced increase with increasing H. 
Remarkably, in spite of the large anisotropy of H C 2, k(H) 
is nearly identical for both H directions at low fields and 
almost saturates at around H s ~ 350 Oe, as shown in the 
inset of Fig. 2(b). Above H s , k(H) becomes anisotropic 
with respect to the field direction and is governed by the 
anisotropy of H C 2- For both H\\ c and ff|| ab, k(H) 
grows gradually with H with a convex curvature, fol- 
lowed by a rapid increase with a concave curvature up to 
H C 2\ there is an inflection point at ~ H c2 /2. 

In superconducting states, the thermal transport is 
governed by delocalized quasiparticles, which shows a 
contrasting field dependence between fully-gapped and 
nodal superconductors, as illustrated in the lower inset 
of Fig. 2(a). In fully-gapped superconductors, the only 
quasiparticles states present at T <C T c are those bound 
to the vortex cores and unable to transport heat, until 
these quasiparticle states within vortices begin to over- 
lap with those within the neighboring vortices. Conse- 
quently, k(H) shows an exponential behavior with very 
slow grows with H in low fields and a rapid increase near 
H C 2 (the dotted line). Such a field dependence is indeed 
reported in fully-gapped superconductors such as Nb [23[ 
and BaNi2As2 [23] • In sharp contrast, the heat transport 
in superconductors with nodes or with a large anisotropy 
in the gap is dominated by contributions from delocalized 
quasiparticles outside vortex cores. The most remark- 
able effect on thermal transport is the Doppler shift of 
the energy of quasiparticles [25], [26| . In the presence of 
line nodes where the density of states has a linear energy 
dependence, k(H) increases in proportion to \[H. Then 
k(H) increases rapidly as soon as the field exceeds the 
lower critical field H c i [26|,|27[ (the solid line). Ifnodeless 
and nodal gaps are mixed in a multiband system without 
interband scatterings, an inflection point emerges in the 
intermediate field regime (the dash-dotted line). 

Now let us discuss the thermal conductivity of LaFePO 
below H s shown in the inset of Fig. 2(b). The nearly 
isotropic k(H) with respect to the field direction defi- 
nitely indicates that the steep increase of k(H) below 
H s is not due to the influence of the first vortex pene- 
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FIG. 2: (color online) (a) Field dependence of k(H) — «(0) 
normalized by k{H c c2 ) - «(0) for H\\ c at T = 0.46 K. Upper 
inset: The same plot in low fields for H || c at 0.46 K (circles), 
1.3 K (triangles) and 1.7 K (squares). Lower inset: Schematic 
field dependence of k(H) — k(0) in superconductors with a full 
gap (dotted line), with line nodes (solid line), and with two 
kinds of gaps with and without nodes (dash-dotted line) . (b) 
Field dependence of [k(H) - k(0)]/[k(H$) - «(0)] for H\\ ab 
plane at 0.46 K. Inset: A comparison of the low-field data for 
if 1 1 ab with those for H || c measured in the zero-field (filled 
circles) and field cooling (FC) conditions (open circles). 



tration field, which is expected to be much smaller than 
H c i ~ 100 Oe, if the demagnetization is taken into ac- 
count 28]. We can also rule out a possibility that the 
steep increase is caused by the remanent Sn-flux inside 
the crystal, because k{T), magnetization, and microwave 
surface impedance measurements [29j show no anomaly 
at T c of Sn (= 3.72 K). Moreover, as shown in the upper 
inset of Fig. 2(a), the low- field steep increase disappears 
at T ~ 1.5 K well below T c of Sn. Thus, we conclude that 
the observed iJ-dependence of k{H) below H s is related 
to the quasiparticle excitations inherent to LaFePO. 

The steep increase and subsequent gradual increase of 
k(H) above H s for both H\\ c and H\\ ab indicate that 
a substantial portion of the quasiparticles is already re- 
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(color online) The same data in the main panel of 

The solid line is a guide 



Fig. 2(a) plotted against sj H/H C 2 
for the eyes. The inset illustrates the extended s-wave (nodal 
s±) gap structure [2(|] in the unfolded Brillouin zone. The 
sign of the gap changes between the solid and dotted lines. 
The small band at the zone corner represents the 3D band 
around the Z point. 



stored at H s , much below H C 2- Such a two-step field 
dependence has been reported in MgB 2 [3(| , NbSe2 [3l[ , 
PrOs4Sbi2 [32| . and URu 2 Si 2 [Hj], providing direct evi- 
dence for the multiband superconductivity. Here H s is 
interpreted as a "virtual upper critical field" that controls 
the field dependence of the smaller gap of the "passive" 
band. Its superconductivity is most likely induced by 
the proximity effect of the "active" bands with primary 
gap. The ratio of the large and smal l gaps is roughly 
estimated to be A^/As ~ ^H^ 2 /H s ~ 6. If we take 
A L - 1.7 k B T c - 13 K, we obtain A s ~ 2 K. This is in 
good correspondence with the fact that the steep increase 
in k(H) below H s is smeared by a similar temperature 
scale (> 1.5 K), as shown in the upper inset of Fig. 2(a). 

The observed multiband superconductivity is compat- 
ible with the band structure of LaFePO, which has been 
examined theoretically and experimentally. The band 
structure calculations show that Fermi surface consists 
of two electronic cylinders centered at the M point and 
two hole cylinders centered at the T point, together with 
a single hole pocket with 3D like dispersion at the Z p oint 
in the Brillouin zone (see the sketch in Fig. 3) [ljj |34| . 
The 3D band is suggested to appear in LaFePO, not in 



LaFeAsO, and has a character of the 3d 3z 2_ r 2 orbital 
which is expected to have a weak coupling to other 2D 
bands [17]. Note that only the 2D cylindrical electron 
and hole bands have been reported by photoemission (3f| 
and dHvA measurements [2ll l34j . 

Nearly isotropic k(H) with respect to the field direc- 
tion below H s shown in the inset of Fig. 2(b) indicates 
that the smaller gap is present most likely in the 3D hole 
pocket. This is consistent with the expected weak cou- 
pling between the 3D and 2D bands. This passive 3D 
band is inferred to be fully gapped because of the fol- 
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lowing reasons. The field dependence of k below H s does 
not show strong \fll dependence expected for line nodes. 
Moreo ver, since th e coherence length of the smaller gap, 
£ s = \J 4>q/(2t:H s ) ~ 100 nm, is comparable to the mean 
free path, it is unlikely that a nodal superconductivity 
can survive against such a "dirty" condition (£ s ~ €}. 

Next we discuss the gap structure of the 2D bands 
from k{H) above H s , where essentially all quasiparticles 
of the 3D band with smaller gap have already contributed 
to the heat transport. As shown in Fig. 3, n(H) for 
H || c increases as ~ VH just above H s to ~ 0.4ii c2 (the 
\/i?-dependence is not clear below H s in our resolution) . 
This ii-dependence and the appearance of the inflection 
point from convex to concave ii-dependence are in sharp 
contrast to the ii-dependence of the simple fully-gapped 
superconductors, in which ii-dependence is always con- 
cave well below T c [23j. Such a convex (sub-linear) ii- 
dependence at low fields appears when the gap is highly 
anisotropic with a large amplitude modulation. On the 
other hand, the concave ii-dependence just below ii C 2 
has never been reported in superconductors with large 
anisotropic gap, such as Tl 2 Ba 2 Cu06+5 [27], CePtsSi 
[11, UPt 3 and LuNi 2 B 2 C Therefore, it is likely 
that at least one of the active 2D bands is fully gapped 
without nodes, although more detailed theoretical inves- 
tigations of the multiband systems are required (39j . In 
fact, the ii-dependence of k(H) with two kinds of gaps 
with and without nodes in the absence of interband scat- 
terings shows an inflection behavior (see the dash-dotted 
line in the lower inset of Fig. 2(a)), which qualitatively 
reproduces the data. This result, along with the finite 



kqq/T observed in our dHvA-available clean crystal, sup- 
ports the nodal superconductivity suggested by the linear 
temperature dependence of the superfluid density [l8[. 
Thus the whole ii-dependence of k above H s implies that 
the 2D bands consist of two kinds; one has nodes and the 
other is fully gapped. 

Finally we discuss the position of the nodes. As can- 
didates of the gap structure with line nodes, the "nodal 
s ± -wave" and "d-wave" symmetries have been proposed 
for LaFePO [2(j. We infer that the cZ-wave can be ex- 
cluded because it possesses line nodes in the 3D band (as 
well as the 2D hole bands), which is unlikely for the rea- 
sons discussed above. The nodal s±-wave structure has 
a nodal gap on the 2D electron band around M point 
and the 2D and 3D hole bands are fully gapped (see the 
sketch in Fig. 3). The gap size of the 2D electron and 
hole bands can be comparable to each other [2(| • There- 
fore, we suggest that the nodal s±-wave structure can be 
the best candidate for the gap symmetry of LaFePO. 

In summary, from the measurements of the thermal 
conductivity, LaFePO is found to be a multigap super- 
conductor with 2D active and 3D passive bands. The 
peculiar field dependence of k provides a stringent con- 
straint on the superconducting gap structure in this sys- 
tem: there exist fully gapped 2D and 3D bands and the 
nodes locate most likely on the other 2D bands. These 
results are consistent with the nodal s±-wave symmetry 
proposed for the superconducting state of LaFePO. 
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